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We study the electronic properties of a twisted trilayer graphene, where two of the layers have 
Bernal stacking and the third one has a relative rotation with respect to the AB-stacked layers. 
Near the Dirac point, the AB-twisted trilayer graphene spectrum shows two parabolic Bernal-like 
bands and a twisted-like Dirac cone. For small twist angles, the parabolic bands present a gap that 
increases for decreasing rotation angle. There is also a shift in the twisted-like Dirac cone with a 
similar angle dependence. We correlate the gap in the trilayer with the shift of the Dirac cone in 
an isolated twisted bilayer, which is due to the loss of electron-hole symmetry caused by sublattice 
mixing in the rotated geometry. Using a tight-binding and a continuum model, we derive an effective 
Hamiltonian which accounts for the relevant low-energy properties of this system. 



I. INTRODUCTION 

The electronic properties of few-layer graphene have 
a remarkable dependence on the stacking arrangement. 
Bilayer graphene with Bernal or AB stacking, which 
is the most energetically favorable structure, has a 
low-energy parabolic spectrum; however, with direct 
or AA stacking, the carriers remain massless with a 
linear dispersion. Rotational disorder in the stacking 
has been experimentally evidenced and theoretically 
studied. Twisted bilayer graphene presents an intriguing 
low-energy behavior as a function of the rotational angle 
between layers, ranging from linear dispersion with a 
Fermi velocity equal to that of monolayer graphene to 
a velocity renormalization for small rotation anglesP^l 
Applying a perpendicular electric field to these bilayer 
phases yields different results: while the A A or twisted 
bilayer graphene remains gapless, for the AB stacking 
there is a gap induced by an external volt age QSHEl 
which is extremely important for device applications. 
Symmetry and changes in the interlayer interaction due 
to the different stacking arrangements are responsible 
for this variety of behaviors. 

Trilayer graphene also presents a remarkable depen- 
dence on the stacking order. Bernal or ABA stacking 
shows an energy dispersion relation that looks like the 
superposition of the monolayer massless linear bands 
and the parabolic massive AB-like bands. Differently 
from the Bernal bilayer, no gap opens in the Bernal 
ABA trilayer under the application of an electric field. 
However, rhombohedral trilayer graphene, with ABC 
stacking, does present a gap with an applied electric 
field.BlffirJl These distinct behaviors can be related to 
the different symmetries of the stackings: the ABA 
trilayer shows mirror symmetry, whereas the ABC 
rhombohedral trilayer does not, having instead spatial 



inversion symmetry, as the Bernal bilayer. Trilayer 
graphene also may present turbostratic disorder, either 
complete, with all layers disoriented, or partial, with 
two layers stacked AB and a third layer showing a 
relative rotation. These twisted trilayers have been 
experimentally found on SiC-grown samples, graphene 
on graphite and in exfoliated few layer grapheneP^HH] 
Some previous theoretical works have partially addressed 
the properties of these twisted trilayer structures, but a 
detailed analysis is still lacking. 

In this work we study a trilayer graphene composed 
of two layers with Bernal stacking and a third one ro- 
tated from an initial ABA position, i.e., an AB-twisted 
layer (ABT). We explore the dependence of the electronic 
properties on the rotation angle 9 employing a tight- 
binding and a continuum model. This allows us to derive 
an effective Hamiltonian which describes the low energy 
behavior of this system. 

Our main results are the following: 

(i) The ABT graphene trilayer shows a combination of 
bilayer Bernal-like parabolic bands with states located 
on the AB-stacked layers, and linear twisted-like bands 
forming a Dirac cone, with states mainly found in the 
twisted layer. 

(ii) The parabolic bands develop a gap that increases 
for diminishing angle. Likewise, a shift of the twisted- 
like Dirac point is observed with a similar angular 
dependence. The velocity of the twisted-like bands is 
renormalized and tends to zero when 9 — > 0. 

(iii) We find that the gap is correlated with a shift of the 
Dirac point of the isolated twisted bilayer. This shift is 
related to the loss of electron-hole symmetry due to the 
mixing of the two sublattices produced by the twisted 
layer. 

(iv) In order to understand the band structure and 
analyze our results, we have compared the outcome of 
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our tight-binding model with that attained by using a 
generalization of the continuum approximation of Lopes 
dos SantosP We obtain that the shift of the Dirac cone 
of the isolated twisted bilayer is erroneously described 
by the continuum model. As a consequence, all the 
low-energy physical properties related to this shift, such 
as the gap and the position of the Dirac point of the 
twisted- like bands in the trilayer, are not appropriately 
described by the continuum model. However, the 
velocity renormalization is suitable described by the 
continuum approximation. 

(v) Equipped with the knowledge gained from the com- 
parison of the tight-binding and the continuum model, 
we derive an effective Hamiltonian which captures the 
low-energy physics of these trilayers, giving a simple 
description of their dispersion relation. 

The paper is organized as follows: we describe the ge- 
ometry of the twisted trilayer and models employed in 
Sec. [II] In Sec. |III| we present our results, which we dis- 
cuss in the light of the continuum model. We finish with 
a summary of our main results in Sec. |IV| 

II. GEOMETRY AND MODEL 

A. Geometry 

We have studied a trilayer graphene composed of two 
layers with Bernal stacking plus a third one with an 
arbitrary relative rotational angle with respect to the 
other two. This twisted layer occupies an outer position 
in the stacking. In order to build a commensurate unit 
cell we have followed a procedure to find coincidence 
lattice points in the crystal similar to that described 
by Campanera et al. for a twisted bilayer P3 First, we 
have three AB-stacked graphene layers (ABA) and then 
we rotate one of the outer layers. We choose a B site, 
i.e., a site where an atom in one layer is exactly at the 
center of the hexagon of the lower layer, as our rotation 
center. We select a site of the layer with coordinates 
r = maj + na.2 and rotate it an angle 8 to an equivalent 
site ti = n&i + ma.2, where ai = (— 1/2, v3/2)ao and 
a.2 = (1/2, v3/2)ao are the graphene lattice vectors; 
n,m are integers, and a = 2.46 A is the graphene 
lattice constant. The unit cell vectors of the trilayer 
twisted cell can be chosen as ti = na\ + ma.2 and 
t2 = — mai + (n + m)a.2- The trilayer unit cell built in 
this way has N = 6(n 2 + ran + m 2 ) atoms and can be 
labeled by the indices (n, ro) , whi ch identify the twisted 
bilayer part of the structure)2il22l The distance between 
layers is set to 3.35 A. 

Figure [l] shows the unit cell for a (2,1) trilayer, where 
the top graphene sheet has been rotated with respect to 
the Bernal-stacked bilayer. This system can be viewed as 
two structures stacked together with a common layer, one 
being a Bernal bilayer graphene and the other a twisted 




FIG. 1. (Color online) Unit cell for the ABT (2,1) trilayer 
graphene. The rotation angle 8 is 28° and the number of 
atoms equals 42. (a) Side view of the unit cell, indicating 
the rotation axis for the twisted layer and a schematic view 
of the hoppings considered in the tight binding model (the 
chosen set of parameters implies hoppings beyond the unit 
cell), (b) Top view of the unit cell. The Bernal-stacked layers 
are shown in red and blue (dark gray) and the twisted later 
in green (light gray). 



bilayer graphene.^ The geometry of the twisted unit cell 
has been extensively studied by Shallcross et al^ and 
MeleP 



B. Tight Binding Model 

We model the trilayer graphene band structure within 
the tight binding approximation including only the p z 
orbitals. Within each layer, we consider a fixed nearest- 
neighbor intralayer hopping parameter 70 = —3.16 eV. 
We take into account the tunneling between adjacent 
layers, neglecting the direct hopping between the outer- 
most onespS For the layer-layer interaction we employ a 
distance-dependent hopping .' 21 " 23 ! Thus, the Hamiltonian 
is given by H = Hi + H 2 + H 3 + H 12 + H23, where H n 
(n = 1, 2, 3) is the Hamiltonian for the individual layer n 
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and H rnn describes the interlayer coupling between con- 
secutive layers to, n: 

H mn = lxe-^- A) c\ Cj + H.c, (1) 

where 71 = —0.39 eV is the nearest-neighbor interlayer 
hopping parameter, d is the interlayer distance, ry is 
the distance between atom i on layer to and atom j on 
the other layer n, and /3 = 3. This value of accurately 
reproduces the dispersion bands ca lculated within a Den- 
sity Functional Theory approach! 21 ! 26 ! 27 ! We let every 
atom in a layer interact with the atoms in the adjacent 
layer located inside a circle of radius GacCi where acc 
is the nearest-neighbor distance between carbon atoms, 
equal to 1.42 A. This takes into account the complexity 
of the unit cell and at the same time breaks the electron- 
hole (e-h) symmetry due to the fact that we are mixing 
the two sublattices. Notice that we employ the same in- 
terlayer Hamiltonian to model the interaction between 
the Bernal-stacked and the twisted layer. 




r < — k — > m r < — k — > m 




FIG. 2. (Color online) Band structures corresponding to four 
instances of ABT trilayer graphene. The corresponding angles 
are indicated in each panel. 



Continuum Model 



III. RESULTS AND DISCUSSION 



We have developed a continuum model for describing 
the twisted trilayer graphene that is a generalization of 
that proposed by Lopes dos Santos^ for twisted bilayer 
graphene. 

The low energy properties of the AB twisted trilayer 
graphene are described in the continuum approximation 
and for angles between 1° and 20°, see Appendix, by the 
following 4x4 Hamiltonian, 
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where v and eo & re the renormalized velocity and the shift 
of the Dirac cone of the twisted bilayer respectively^ and 
g±=g(k± ^jp) with g(k)=k x — ik y , AK being the vector 
connecting the Dirac cones of the rotated layers. The 
basis for this Hamiltonian arc the non-dimcr atoms A\ 
and B 2 of layers 1 and 2 respectively, and the two basis 
atoms A3 and B3 of layer 3. From the diagonalization 
of the previous Hamiltonian we obtain the following four 
bands, 
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Notice that the continuum model describes the bands 
corresponding to a given valley of the original graphene 
layer. Interchanging g + by g_ in Eq. [2] yields the results 
for the other valley, for which the position in k space of 
the parabola and the Dirac cone are swapped. In order 
to compare to the TB results, the bands for both valleys 
should be superposed. This is equivalent to substitute 
g± by g in Eq. ||| 



Figure [2] shows the tight-binding band structures for 
four instances of ABT twisted trilayers with angles com- 
prised between 1.35° and 5.09°. Two linear bands corre- 
sponding to the twisted layers and the parabolic bands 
stemming from the Bernal-stacked layers are clearly seen 
for all the cases depicted. For the largest angle shown in 
Fig. [2] (a), 5.09 °, the four bands are practically almost 
degenerate at zero energy and in the scale of energies of 
the plot the spectrum looks like e-h symmetric. How- 
ever, for smaller angles, a gap between the two parabolic 
bands and a shift of the Dirac point in the linear bands 
can be discerned. In all the depicted twisted angles, the 
lowest energy parabolic band is pinned at zero energy. 
This band is very weakly coupled and has a non-bonding 
character; for this reason its energy at the Dirac point 
does not vary with the rotation angle. The Dirac point 
and the minimum of the other parabolic band have en- 
ergy shifts that do depend on the twist angle. Besides, 
we find that the degree of e-h asymmetry increases dra- 
matically with diminishing angle, see Fig. [2] panels (b) 
to (d). 

Velocity renormalization. The Dirac-like bands clearly 
show a Fermi velocity reduction when the twisted angle 
decreases. This is patent for the smaller angle depicted, 
1.35°, for which a reduction of near 90% is observed 
(Figj2](d)). In order to quantify the velocity renormal- 
ization, in Fig. [3] we compare the velocities of the twisted 
Dirac-like bands in the trilayer graphene with those of the 
twisted bilayer case, both obtained in the tight-binding 
approximation. We see that the behavior of the two ve- 
locities is similar up to 1°. In both cases the band velocity 
is strongly renormalized for rotation angles below 5°, and 
tends to zero when 8 — > 0. 

Within our model TB Hamiltonian we find for the 
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FIG. 3. (Color online) Fermi velocities of the trilayer graphene 
twisted-like bands and the twisted bilayer graphene normal- 
ized to the Fermi velocity of monolayer graphene vf calcu- 
lated within the TB approximation. In the case of the tri- 
layer, we stop the plot of the velocity at the angle for which 
a gap develops at the Dirac point. Inset: Fermi velocities 
of the trilayer twisted-like bands and the isolated twisted bi- 
layer calculated with the continuum model. Several zeros of 
the velocities can be distinguished. 



twisted bilayer that the velocity vanishes at an angle near 
0.8°, and then the velocity increases again for smaller an- 
gles. It has been proposed that for even smaller angles 
new zeros in the velocity should appearP As in TB cal- 
culations smaller twisted angles correspond to larger unit 
cells, we do not find a second zero in the Fermi velocity 
because of the numerical limitations to the angles we can 
study. In the twisted trilayer case the situation is differ- 
ent: we find that for angles smaller than 1.16°, where the 
velocity is near zero, the conduction and valence bands 
couple and a small anticrossing gap of the order of 10~ 4 7o 
appears. As we will show below, this anticrossing occurs 
because at small angles the layers are strongly coupled, 
resulting in wave functions with a mixed weight in the 
three layers. By restricting the tunneling between the 
twisted bilayer and the extra layer only to vertical hop- 
pings, we reduce the mixing between layers and this an- 
ticrossing disappears. 

In order to analyze our results we have computed the 
band structure of the trilayer twisted bilayer in a contin- 
uous model. In the inset of Fig. [3] we plot the velocity of 
the Dirac point obtained with the continuous model for 
the twisted bilayer and trilayer graphene. As in the tight- 
binding results, we see that the velocity renormalization 
in the bilayer and trilayer cases is practically the same. 
The continuum approach allows for the exploration of 
smaller angles than the TB model. Thus, it is found that 
the zeros in the velocity appear at the same angle, both 
in the bilayer and in the trilayer cases (see inset of Fig 3| , 
and they correspond to quasi-confined states in the AA- 



stacked regions of the twisted bilayer .EHSI Contrary to the 
TB results, there is no gap in the Dirac-like bands at 
small rotation angles in the continuum approximation. 
As we commented above, this is because within this ap- 
proach the tunneling between layers 1 and 2 is treated 
as vertical, and this yields a zero gap also in the TB 
calculation. 
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FIG. 4. (Color online) Electronic probability densities corre- 
sponding to the (a) parabolic and (b) linear bands distributed 
in the three layers. Red, blue and green colors and numbers 
1, 2, 3 code the three layers from bottom to top, namely (1) 
outer Bernal layer, (2) middle Bernal layer, and (3) twisted 
layer, as depicted in panel (c). Letters a, b, c label the shifted 
parabolic, linear, and zero-energy pinned parabolic bands re- 
spectively, as indicated in the schematic band structure shown 
in panel (d). 



Layer distribution. The spatial distribution of the car- 
riers corresponding to different bands gives an idea of 
the degree of coupling in the system. The carriers cor- 
responding to the band pinned at zero energy, labeled 
c in Fig. Qd), are situated in the outer Bernal layer 
(1), independently of the rotation angle. The behavior 
of the other bands is more complex. In panels (a) and 
(b) of Fig. [i] we present the electronic probability den- 
sities of the shifted parabolic band a and of the linear 
twisted-like bands b, distributed in the three layers, see 
Fig. |4Jd). We have plotted the probability densities of 
the eigenstates close to the Dirac point, checking that the 
two states corresponding to the Dirac cone have the same 
distribution. The shifted parabolic band a corresponds 
to states mainly located in the inner Bernal layer, labeled 
2, especially for large rotation angles. The probability of 
being in the twisted layer (3) increases with diminishing 
angle, most notably below 5°. Likewise, the carriers cor- 
responding to the linear twisted-like band, band b in Fig. 
|4jd), are mostly located in the twisted layer (3), but the 
probability of being in the other two layers increases be- 
low 5°, being almost equally distributed among the three 
layers for angles smaller than 1°. Thus, for smaller an- 



5 



gles we obtain that the carriers corresponding to bands 
a and b have an increasingly mixed spatial distribution, 
evidencing an enhancement of the coupling. This angle 
of 5° is where the renormalization of the velocity for the 
twisted trilayer carriers is noticeable; as it is well known, 
this is also a sign of the increased coupling, that we also 
relate here to the spatial distribution of the electrons cor- 
responding to the low-energy bands. 
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FIG. 5. (Color online) Shift of the Dirac cone (triangles) and 
gap between the parabolic bands (circles) for ABT trilayer 
graphene along with the shift for the Dirac cone for the iso- 
lated twisted bilayer (crosses) calculated with the TB model. 
Inset: Shift of the Dirac cone for the ABT trilayer obtained 
with the continuum model. 

Energy gap and shift. In Fig. [5] we show the shift 
of the Dirac cone and the gap between the Bernal-likc 
parabolic bands found for the ABT trilayer calculated 
with the TB model. For comparison, we also plot the 
shift for the Dirac point in an isolated twisted bilayer 
graphene calculated with our tight-binding model. The 
gap increases with diminishing angle, following closely 
the shift of the Dirac cone for the isolated bilayer. 

In the continuum approximation, there is also a gap 
in the parabolic bands that is equal to the shift of the 
Dirac cone eb, Eq. [3j However, the dependence of this 
gap with the rotation angle is at variance with the tight- 
binding results. In the inset of Fig. [5] we plot eo as 
function of the rotation angle. The shift is obtained by 
diagonalizing the full continuum Hamiltonian in a large 
enough plane wave basis, see Appendix. In the contin- 
uum approximation the gap has an almost constant value 
equal to jf I (27r-\/37o) for angles in the range 4° and 20°, 1 
and decreases to zero at smaller angles. The oscillations 
observed in the continuum results are related to those 
observed in the renormalized velocity (inset of Fig. [jjj. 
We attribute the discrepancy between the tight-binding 
and continuum results to the fact that the continuum 
model can only cope with the layer-layer interactions lo- 
cally. The vertical tunneling neglects the effect of odd- 



numbered carbon atom rings among layers considered in 
a TB model with hopping beyond nearest neighbors. As 
a consequence, the continuum model does not appropri- 
ately describe the electron-hole asymmetry of the system. 

Fig. [5] also shows that around 5°, the two parabolic 
bands cross, so for large angles the non-bonding band 
pinned at zero has a higher energy than the shifted 
parabolic band at the Dirac point. This change in the 
position of the two parabolic bands coincides with the 
change in the spatial distribution of the states corre- 
sponding to the linear twisted-like band and the shifted 
parabolic band. 

Continuum effective model. As the continuum model 
outlined in the previous Section is not valid for the de- 
scription of the low-energy spectrum of twisted trilayer 
graphene, we modify it accordingly, deriving an effec- 
tive Hamiltonian for the band structure near the Dirac 
point. This Hamiltonian describes the parabolic Bernal- 
like bands and the linear Dirac cone corresponding to the 
twisted layer: 
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is the wavevector measured from the Dirac point, and 
eo, ei are small onsite energy terms that break the e-h 
symmetry, giving rise to the shift of the twisted cone and 
the opening of the gap respectively. This model yields a 
decoupled linear Dirac cone stemming from the twisted 
layer and the parabolic bands from the Bernal stacking. 

The difference between this Hamiltonian and that 
given by Eq. [2] amounts to the diagonal terms eo an d 
ex, which are not equal in value and angle-dependent. 
As discussed above, this is due to the loss of e-h sym- 
metry, not described appropriately by the continuum ap- 
proximation. Therefore, we consider them as parameters 
given by the TB calculation. 

Finally, the effect of an external applied bias V be- 
tween the outer layers 1 and 3 is described by adding 
de corresponding diagonal terms to the effective Hamil- 
tonian, Eq. [4] This does not open a gap in the Dirac 
cone, but produces a shift in the Dirac point and modi- 
fies the gap between the parabolic bands. As an example, 
adding — V/2 to the first diagonal element and +V/2 to 
the diagonal of the 2x2 box corresponding to the twisted 
layer, induces a shift V/2 in the Dirac cone and modifies 
the parabolic gap to ei — V/2. Thus, the application of 
a bias may close the gap of the parabolic bands and, for 
negative V, it can shift down the Dirac cone. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have shown that the low-energy spec- 
trum of the twisted trilayer graphene composed by two 
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Bernal-stacked layers and a third rotated graphene sheet 
is a combination of the bilayer Bernal-like parabolic 
bands and linear Dirac cone reminiscent of the twisted 
bands. The parabolic bands show a gap that gets larger 
for smaller angles. There is also a shift in the Dirac cone 
that has a similar angular dependence. The velocity of 
the twisted-like bands is renormalized for small angles, 
tending to zero when the rotation angle tends to zero. 
We find that the gap is correlated with the shift of the 
Dirac cone in an isolated twisted bilayer. The shift is 
due to the loss of the electron-hole symmetry produced 
by the mixing of the two sublattices. 

One of the parabolic bands has a non-bonding char- 
acter, pinned at zero energy at the Dirac point, and its 
carriers are located exclusively in the outer Bernal layer. 
The other bands have a carrier distribution that depends 
on the rotation angle. For larger angles, the parabolic 
band is mostly located on the inner Bernal layer, whereas 
the carriers of the linear bands are mainly on the twisted 
layer. However, for angles below 5° there is an increasing 
weight on the twisted layer for the parabolic band, and 
the linear bands, initially located on the twisted layer, 
gain weight in the other two AB-stacked layers. 

The change in localization indicates an increase in the 
effective coupling between layers, and coincides with a 
change in the relative position of the two parabolic bands, 
i. e., a zero gap point. This also concurs with the on- 
set of the velocity renormalization, consistent with the 
aforementioned enhancement of the effective interlayer 
interaction. 

We have compared our results obtained with a tight- 
binding model with those derived from a continuum ap- 
proximation. While the velocity renormalization it is cor- 
rectly given by the continuum model, the shift of the 
Dirac cone is not well described by this approximation, 
neither with respect to the sign nor to the angular de- 
pendence. We thus derive an effective model for the low 
energy physics of the system, allowing for a simpler de- 
scription of its properties. 



B2 atom (Fig. [T|). Using a superlattice unit cell de- 
fined by the basis vectors ti = (i + l)ai + i&2 and 
t 2 = — ia.i + (2i + l)a 2 , where i is a positive integer, 
we get a commensurate rotation angle defined by 
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We are interested in the low energy spectrum near 
a Dirac point, K=^(1,0), in layer 1 and 2. Under 
the rotation, the equivalent Dirac point in layer 3 is 
K e = ^ (cos 6, sin 6). Setting AK=K e - K, we can 
define wave functions with the same wavevector k in the 
three layers that refer to the same plane wave states in 
the original lattice: the Dirac points occur at — AK/2 in 
layers 1 and 2 and at AK/2 in layer 3. With this the 
intralayer part of the Hamiltonian takes the form 

H cl +H c2 + H c3 =v F J2 4,^^ ■ (k + AK/2) c hK 
+v F £ c 2,Ka (Ta ^ ' (k + AK/2)c 2W 
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with cr e — e t0a *> 2 (a x ,(Ty)e l9cr ^/ 2 . In this expression the 
operator c^~ k creates and electron in layer i on atom a 

and momentum k and i>i?=\/3/27o&o- 

The tunneling between layers 1 and 2, which are AB- 
stacked, is described by the interlayer Hamiltonian 
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where the hopping matrix for the AB stacking is 
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Appendix A: Continuum model for the twisted 
trilayer graphene 



We describe the coupling between the twisted layers 2 
and 3 with the approximation introduced by Lopes dos 
SantosfQ We consider a periodically modulated local tun- 
neling 
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where the vector G is summed over the reciprocal lattice 
vectors of the supercell. In the continuum approxima- 
tion, only the largest t^{G) are retained in the hopping 
Hamiltonian. They correspond to the reciprocal lattice 
vectors 0, Gi and — Gi — G 2 and have the form 



We consider a twisted trilayer graphene formed by 
three graphene monolayers, with two of them (1 and 2) 
A_B-stacked, and the third one rotated a commensurate 
angle 9 about a perpendicular axis passing through a 
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FIG. 6. (Color online) Band structure of twisted trilayer 
graphene as obtained in the continuum approximation. In 
the calculation 71=0.1270. 



with / = e l2jr / 3 and 7=71/3. In the previous expressions 
Gi and G2 are the primitive vectors of the reciprocal 
space defined by the supercell lattice vectors ti and t2- 

The band structure is obtained by diagonalizing the 
Hamiltonian 
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in a plane wave basis. The number of plane waves needed 
to obtain an accurate spectrum increases with the size 
of the supercell, and therefore the calculation is more 
demanding numerically for smaller twisted angles. We 
have checked the convergence of our results with respect 



to the number of plane waves used. In Fig. [6] we plot 
the low energy bands, as obtained by diagonalizing the 
continuum Hamiltonian for three values of the twisted 
angle. The spectrum consists of a Dirac cone located 
mainly on the twisted layer and two gapped parabolic 
bands with states in the AS-stacked bilayer. The con- 
tinuum results are the following: (i) the velocity of the 
Dirac cone v and the curvature of the parabolic bands 
decrease with diminishing rotation angle; (ii) there is an 
energy shift of the Dirac cone, eo, equal to the gap be- 
tween the parabolic bands; (iii) the top of the parabolic 
valence band is pinned at zero energy; (iv) the shift and 
gap decrease when the twisted angle decreases. We have 
checked that eo and v coincide with the energy shift and 
the velocity rcnormalization occurring in twisted bilayer 
graphene From the previous observations we can define 
a Hamiltonian which describes the low energy physics of 
twisted trilayer graphene for angles between 1° and 20°, 
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where g±=g(k.± =Sp) and g(k)=k x — ik y . This Hamilto- 
nian can be further simplified by eliminating the dimcr 
states in the AB stacked layers,^ 
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